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Abstract 

We develop a very general theory on the regularized approximate invisibility 
cloaking for the wave scattering governed by the Helmholtz equation in any space 
dimensions N >2 via the approach of transformation optics. There are four major 
ingredients in our proposed theory: 1). The non-singular cloaking medium is ob- 
tained by the push- forwarding construction through a transformation which blows 
up a subset in the virtual space, where £: <C 1 is an asymptotic regularization 
parameter. will degenerate to Kq as £ ^ +0, and in our theory Kq could be any 
convex compact set in M^, or any set whose boundary consists of Lipschitz hyper- 
surfaces, or a finite combination of those sets. 2). A general lossy layer with the 
material parameters satisfying certain compatibility integral conditions is employed 
right between the cloaked and cloaking regions. 3). The contents being cloaked could 
also be extremely general, possibly including, at the same time, generic mediums 
and, sound-soft, sound-hard and impedance-type obstacles, as well as some sources 
or sinks. 4). In order to achieve a cloaking device of compact size, particularly for 
the case when is not "uniformly small" , an assembly- by-components, the (ABC) 
geometry is developed for both the virtual and physical spaces and the blow-up 
construction is based on concatenating different components. 

Within the proposed framework, we show that the scattered wave field Ue cor- 
responding to a cloaking problem will converge to ixo as e ^ +0, with uq being 
the scattered wave field corresponding to a sound-hard Kq. The convergence result 
is used to theoretically justify the approximate full and partial invisibility cloaks, 
depending on the geometry of Kq. On the other hand, the convergence results are 
conducted in a much more general setting than what is needed for the invisibility 
cloaking, so they are of significant mathematical interest for their own sake. As for 
applications, we construct three types of full and partial cloaks. Some numerical 
experiments are also conducted to illustrate our theoretical results. 
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1 Introduction 

This paper is concerned with the invisibility cloaking for the wave scattering governed 
by the Helmholtz equation via the approach of transformation optics [161 El [36], 
which is a rapidly growing scientific field with many potential applications. We refer to 
[HI [m dni [331 ISHl [iQj and the references therein for discussions of the recent progress on 
both the theory and experiments. 

Let ft and D be two bounded Lipschitz domains in R^, N > 2, such that D ^ ft. 
Let cr = cr(x) = (a^-^(x)) G R^^^, x E R^, be a symmetric-matrix valued measurable 
function such that, for some A, < A < 1, we have 

A||Cf < ^'^'(^)^^^J ^ for any ^ G R^ and for a.e. x G R^. (1.1) 



Let q — qi + iq2 — q{x), x G R , be a complex- valued bounded measurable function 
with real and imaginary parts qi and q2 respectively, such that, for some A, < A < 1, 
we have 

qi{x) > A, q2{x) > for a.e. x G R^. (1.2) 

Furthermore, we assume that q{x) = go 1 and a^^ (x) = := 6ij for x G R^\fi, 
where 6ij denotes the Kronecker delta function. In the following, (1.1) and (1.2) will be 



referred to as the regular conditions on a and q, and A is called the regular constant 

Next, we introduce the time-harmonic wave scattering governed by the Helmholtz 
equation whose weak solution is u = u{x, d, fc), x G R^, where d G k G R+, 

(1.3) 



div(crVi/) + k^qu = in 
u{x, d, k) — e^^^'^ satisfies the radiation condition. 



The last statement in ( |1.3[ ) means that if one lets u^{x, d, k) = u{x^ d, k) — e^^^'^, then 
lim f du{x) | ^ q ^ ^ n n 

r^oo \ or J 



In the physical situation, (1.3) can be used to describe the time-harmonic acoustic 
scattering due to an inhomogeneous acoustical medium (fi; a, g) located in an otherwise 
uniformly homogeneous space (R^\f^; ctq, g'o)- ^ and respectively, denote the den- 
sity tensor and acoustic modulus of the acoustical medium, and u{x^ denotes the wave 
pressure with U{x,t) :— u{x)e~'^^^ representing the wave field satisfying the scalar wave 
equation 

q{x)Uu{x,t) -J2q^. (^''^^^a^^^^'^) ) ^ ° in R^ X R. 
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The function u'^{x) := e is an incident plane wave with k denoting the wave number 
and d E denoting the impinging direction. u{x) is cahed the total wave field and 

u^{x) is called scattered wave field, which is the perturbation of the incident plane wave 
caused by the presence of the inhomogeneity (fi; a, q) in the whole space. Indeed, it is 
easily seen that if there is no presence of the inhomogeneity, will be vanishing. For the 
particular case with N — 2^ ( |1.3[ ) can also be used to describe the transverse-electric (TE) 
polarized electromagnetic (EM) wave propagation with the presence of an infinitely long 
cylindrical EM inhomogeneity (fi;cr, (see, e.g., [7J). In this case, e \— 5ftg, 7 := k^q 
and fi '.— denote, respectively, the electric permittivity, conductivity and magnetic 
permeability, where 5ft and denote taking the respective real and imaginary parts. We 
refer to j9l[32] for related physical background. In the rest of the paper, in order to ease 
the exposition, we stick to the terminologies with the acoustic scattering. 



We recall that by a weak solution to (1.3) we mean that u G i7/^^(R^) and that it 
satisfies 



/ 



aVu • - k^quip = for any ip E C^(M 



The limit in (1.4) has to hold uniformly for every direction x — x/\\x\\ G E>^~^ and 
is also known as the Sommerfeld radiation condition which characterizes the radiating 
nature of the scattered wave field (cf. |^). There exists a unique weak solution 
u{x^d^k) = u~x^ + '^^XRA^\n ^ ^/oc(^^) ( |1-3| ), and we refer to Appendix in [27] 
for a convenient proof. We remark that, if the coefficients are regular enough, (1.3) 
corresponds to the following classical transmission problem 



N 



\ a^^ ^—u (x,d,k)]^k^qu (x, d. A:) = 
dxi V oxj J 

Ai/+(x, d, k) + k^u^{x, d,k) = 



N 



; du 



U (x) = U^{x), ^ i^i^^^ -Q )(^) = • V'^^)(5 

ij=l ^ 

i/+(x, d, k) = e'^^'"^ + u%x, d, k) 



lim r 2 f du{x) ^ 
r^cx) \ or 







X G M^\n, 

X G dVL, 

X G M^\n, 

r — llxIL 



(1.5) 



where v — is the outward unit normal vector to dO.. 



Furthermore, u{x) admits the following asymptotic development as ||x|| +oc 



u{x^ d, k) 



^ikx-d 



ik\\x\\ 



N-l 

\x\\ 2 



-, d, 



+ 



2 



(1.6) 



In (1.6), Uoo{x,d,k) with x := ^ S is known as the far-field pattern or the 

scattering amplitude^ which depends on the impinging direction d and wave number k of 
the incident wave u'^{x) := e^^^'^, observation direction x, and obviously, also the under- 
lying scattering object (fi; a, q). In the following, we shall also write Uoo(x^ d; (fi; a, q)) to 
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indicate such dependences, noting that we consider k to be fixed and we drop the depen- 
dence on k. An important inverse scattering problem arising from practical applications 
is to recover the medium (fi;cr, by knowing of Uoo{x,d). This inverse problem is of 
fundamental importance to many areas of science and technology, such as radar and 
sonar, geophysical exploration, non-destructive testing, and medical imaging to name 
just a few; see [91 [TSJ and the references therein. In this work, we shall be mainly con- 
cerned with the invisibility cloaking for the inverse scattering problem, which could be 
generally introduced as follows. 

Definition 1.1. Let and D be bounded Lipschitz domains such that D ft\D 
and D represent, respectively, the cloaking region and the cloaked region. Let F and 
be two subsets of E>^~'^. {^\D; ac, Qc) is said to be an (ideal/perfect) invisibility cloaking 
device for the region D if 

Uoo{x,d;{Q;ae,qe)) = for x G F, d G F^, (L7) 
where the extended object 



acQc in Q\D, 



with (D^aa^Qa) denoting a target medium. If F = F' = S^"-*^, then it is called a full 
cloak^ otherwise it is called a partial cloak with limited apertures F of observation angles, 
and F' of impinging angles. 



By Definition 1.1, we have that the cloaking layer {Vt\D]ac^qc) makes the target 
medium (D; a^, qa) invisible to the exterior scattering measurements when the detecting 
waves come from the aperture F' and the observations are made in the aperture F. 

One efficient way of constructing the invisibility cloak that has received significant 
attentions in recent years is the so-called transformation optics [IGj [TTl |22l [36]. By 
taking advantage of the push-forward properties of the material parameters a and 
the transformation optics approach via a blow-up transformation in constructing an 
(ideal) invisibility cloak can be simply described as follows. Let (fi;(Jo,g'o) be selected 
for constructing the cloaking device, and let P G be a point. ctq, go) lives in the 

so-called virtual space. Suppose that there exists a transformation F which blows up the 
point P to an open subset D within fi. The homogeneous virtual space is then pushed- 
forward to form the cloaking layer (Q\D;ac,qc)- The cloaking layer together with a 
filling-in target medium {D;aa,qa) forms the cloaking device, which lives in the so-called 
physical space. Due to the transformation invariance of the Helmholtz equation, it can 
be heuristically argued that the scattering amplitude in the physical space is the same as 
the scattering amplitude in the virtual space. Since the scatterer in the virtual space is 
a singular point P, whose scattering effect is negligible, this implies that the scattering 
amplitude in the physical space is also vanishing. Here, we would like to emphasize 
that from a practical viewpoint, the target medium should be arbitrary or as general 
as possible, and this viewpoint shall be adopted throughout our current study. The 
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blow-up-a-point construction yields singular cloaking materials, namely, the material 
parameters violate the regular conditions. The singular media present a great challenge 
for both theoretical analysis and practical fabrications (cf. [HI [29]). In order to avoid the 
singular structure, several regularized constructions have been developed. In [lOl |T2l ES] , 
a truncation of singularities has been introduced. In [JOlEllES], the 'blow-up-a-point' 
transformation in [l71|22l[36j has been regularized to become the 'blow-up-a-small-region' 
transformation. In the current study, we shall adopt the latter one for the construction 
of our cloaking device. Nevertheless, as pointed out in [19j, the truncation-of-singularity 
construction and the blow-up-a-small-region construction are equivalent to each other. 
Hence, all the obtained results in this work equally hold for the truncation-of-singularity 
construction. Instead of ideal/perfect invisibility, one would consider approximate/near 
invisibility for a regularized construction; that is, one intends to make the corresponding 
scattering amplitude due to a regularized cloaking device as small as possible depending 
on an asymptotically small regularization parameter e G M+. This is the main subject 
of study for the present paper. 

Due to its practical importance, the approximate cloaking has recently been ex- 
tensively studied. In [5l [21], approximate cloaking schemes were developed for EIT 
(electric impedance tomography) which might be regarded as optics at zero frequency. 
In [6l [U [201 [231 [2H1 ISSl [HI [35] , various near-cloaking schemes were presented for scalar 
waves governed by the Helmholtz equation. In all the aforementioned work, the con- 
structions of the cloaking layer {^\D] a^, q^) are based on blowing up a uniformly small 
neighborhood of a singular point P; namely, degenerates to the single point P 
as 5 ^ +0. In order to stabilize and enhance the accuracy of approximation of the 
near-cloaks, various mechanisms have been developed in those literatures. Particularly, 
we would like to note that, in [20], it is shown that the regularized approximate cloak is 
unstable due to the existence of cloak-busting inclusions, and the authors propose to in- 
corporate a special lossy layer to stabilize the approximation. A different lossy layer was 
proposed and investigated in [231 [28]. The cloaking of impenetrable obstacles, which 
could be taken as lossy mediums with extreme material parameters, were considered 
in [3 [SI [1] and [26]. Also, we would like to point out that, in all those studies, the 
approximate full invisibility cloaks were obtained. 

In the present work, we develop a very general theory on the regularized approximate 
invisibility cloaking for the wave scattering governed by the Helmholtz equation in any 
space dimensions N > 2 via the approach of transformation optics. First, the non- 
singular cloaking medium is obtained by the push-forwarding construction through a 
transformation which blows up a subset in the virtual space, with degenerating 
to as 5 ^ +0. In our theory, could be very general. It could be any convex 
compact subset in M^, or any set whose boundary consists of Lipschitz hypersurfaces, 
or a finite combination of those sets. For example, in M^, it could be a single point, 
or a line segment, or a bounded planar subset. This includes all the existing studies 
in the literature by blowing up 'point-like' regions as a very special case. Second, in 
order to stabilize the approximation process, a lossy layer with the material parameters 
satisfying certain mild compatibility integral conditions is employed right between the 
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cloaked and cloaking regions. The lossy layer is also very general and could be variable 
and even be anisotropic. Third, the proposed cloaking scheme is shown to be capable 
of nearly cloaking an very general content, possibly including, at the same time, generic 
passive mediums, and sound-soft, sound-hard, and impedance-type obstacles, and some 
active sources or sinks as well. Finally, in order to achieve a cloaking device of compact 
size, particularly for the case when is not 'point-like', assembled- by-components 
(ABC) geometry is developed for both the virtual and physical spaces and the blow- 
up construction is based on concatenating different components. Within the proposed 
framework, we show that the scattered wave field corresponding to a cloaking problem 
will converge to uq SiS £ ^ +0, with uq being the scattered wave field corresponding to a 
sound-hard Kq. The convergence result is used to theoretically justify the approximate 
full and partial invisibility cloaks, depending on the geometry of Kq. On the other 
hand, the convergence results are conducted in a much more general setting than what 
is needed for the invisibility cloaking, so they are of significant mathematical interest for 
their own sake. As for applications, we construct three types of full and partial cloaks. 
Some numerical experiments are also conducted to illustrate our theoretical results. 

It is interesting to note that in addition to the blow-up-a-single-point construction, 
the cloaking constructions based on blowing up an arc curve or a planar rectangle are also 
proposed and investigated in [IHlEl], and they respectively yield the so-called electro- 
magnetic wormholes and carpet-cloaking. As discussed earlier, the regularized blow-up- 
a-single-point construction, namely the blow-up-a-small-region construction, has been 
extensively studied in the literature. Using the general framework developed in the 
present work, one can easily construct the regularized electromagnetic wormholes and 
carpet-cloaking by employing non-singular materials. In this paper, we focus entirely on 
the transformation optics approach in achieving the cloaks. We would like to mention 
in passing other promising cloaking techniques which we did not consider in the present 
study including the one based on anomalous localized resonance [31 [3T], and another one 
based on special (object-dependent) coatings [2]. 

The rest of the paper is organized as follows. In the next section, we present the 
general blow-up construction of the proposed regularized cloaks and give some relevant 
discussions. Section 3 is devoted to the convergence analysis in the virtual space. Section 
4 is on the application of the results obtained in Section 3 to the construction of full 
and partial cloaks in the physical space. In Section 5, we develop the ABC-geometry for 
both the virtual and physical spaces, and construct three types of full and partial cloaks 
that are new to the literature. Finally, in Section 6, we give some numerical simulations. 

2 General construction of the regularized cloaks 

In this section, we shall give the general construction of a regularized cloaking device 
via the transformation optics approach based on a blow-up mapping between the virtual 
and the physical spaces. The main purpose of this section is to pave the way for our 
convergence analysis study in the virtual space that shall be conducted in the next 
section. We first give a definition of an admissible acoustic configuration. 
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For any x G M^, > 2, we denote x = {x'^xn) G M^~-'^xR and x = (x^^, xat-i, xn) G 
R^~^ X M X R. For any r > and any x E M^, denotes the Euclidean ball con- 

tained in with radius r and center x, whereas S^(x') denotes the Euclideean ball 
contained in R^"-*^ with radius r and center x^ Moreover, = Br{0) and = S^(0). 
Finally, for any £; C R^, we denote Br(E) = |Jxg£; ^r(^). 

Definition 2.1. We say that K C R^ is a scatterer if K is compact and G = R^\i^ is 
connected. 

A scatterer K C Br, for some > 0, is regular if the immersion ^ 
L'^{Bji^i\K) is compact. 

We say that a scatterer K is Lipschitz-regular if, for some positive constants r, Li 
and L2, for any x G there exists a bi-Lipschitz function • Br{x) R^ such that 
the following properties hold. First, for any zi, Z2 ^ Br{x) we have 

Ll\\zi - Z2\\ < \\^x{zi) - ^x{z2)\\ < L2\\zi - Z2\\. 

Second, = and ^x{dK D Br{x)) C tt = G R^ : Vn = 0}- 

A scatterer K is said to be Lipschitz if, for some positive constants r and L, the 

following assumptions hold. 

For any x G dK, there exists a function (f : R^"-*^ — > R, such that ip{0) = and 

which is Lipschitz with Lipschitz constant bounded by L, such that, up to a rigid change 

of coordinates, we have x = and 

Br{x) ndKc{ye Br{x) : vn = ^{y)}. 

We say that x G dK belongs to the interior of dK if there exists 5, < 5 < r, 
such that Bs{x) H dK = {y ^ Bs{x) : i/n = ^{y')}- Otherwise we say that x belongs 
to the boundary of dK. We remark that the boundary of dK might be empty and 
that, if X G dK belongs to the interior of dK, then K may lie at most on one side of 
dK, that is Bs{x) DK = Bs{x) D dK, or Bs{x) n K = {y e Bs{x) : yN > ^{y')}, or 
Bs{x) DK^iye Bs{x) : yN < ^{y^)}- 

For any x belonging to the boundary of dK, we assume that there exists another 
function (fi : R^~^ R, such that (pi{0) = and which is Lipschitz with Lipschitz 
constant bounded by L, such that, up to the previous rigid change of coordinates, we 
have X = and 

Br{x) ndK = {ye Br{x) : yN = ^{y'), yN-i < Mv'')}- 

Finally, for any x G dK, let ei(x), . . . , eAr(x) be the unit vectors representing the 
orthonormal base of the coordinate system for which the previous representations hold. 
Then we assume that eAr(x) is a Lipschitz function of x G dK, with Lipschitz constant 
bounded by L, and eN-i(x) is a Lipschitz function of x, as x varies on the boundary of 
dK, with Lipschitz constant bounded by L. 

Properties of Lipschitz scatterers are thoroughly investigated in [30l Section 4]. Let 
us just notice that a Lipschitz scatterer K is Lipschitz-regular. Furthermore, a Lipschitz- 
regular scatterer K C Br is regular and the immersion W^''^{Br^i\K) L'^{dK) is 
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compact. Notice that for a connected component K of K with empty interior, with a 
shght abuse of notation, with dK we denote two copies of {K^, and by L?{dK) 
we denote the couple {u^^u~) G L?{K^) x L?{K~). In such a way we can define the 
trace of a function u G W^''^{B^j^i\K) on both sides of K. 

Definition 2.2. We fix > 0. An acoustic configuration C = (Ki, ^3, <§, a, /i, H) 

is admissible if the fohowing assumptions hold. 

There exist three scatterers i^i, K2, which are pairwise disjoint (possibly some or 
all of them may be the empty set) and such that K2 is regular and Ks is Lipschitz-regular. 
We set K = K1UK2U Ks. 

Let 5 = 5i + Z52 = s{x), X G dK^, be a complex- valued bounded T/^"-*^ -measurable 
function, with real and imaginary part si and 52 respectively, such that 

si{x) < and S2{x) > for n^'^-di.e. x G 8X3. 

Again, on a connected component K of K3 with empty interior, s G L^(dK) means 
(5+, 5-) G L^{K^) X L^(K-). 

Let a = cr(x), x G M^\i^, be an x symmetric matrix whose entries are real- 
valued measurable functions such that, for some A, < A < 1, we have 

A||Cf < (r{x)^ • C < X~^Uf for any ^ G and for a.e. x G R^\K. 

Let q — qi + iq2 — q{x), x G R^\K, be a complex- valued bounded measurable 
function, with real and imaginary part qi and q2 respectively, such that, for some A, 
< A < 1, we have 

qi{x) > A and q2{x) > for a.e. x G R^\K. 

Furthermore, we assume that 

a(x) = I and q{x) = 1 for any x outside a compact subset. 

The source term {h, H) G L^(R^\K, C x C^) and it has compact support. 
Finally, we require that the following scattering problem has only a trivial weak 
solution 

' diY{aVu)^Pqu = mR^\K, 
u = on dKi , 

aVu-u^O ondK2, (2.1) 

(jVu ' v — su — {) on 5^3, 

u satisfies the radiation condition, 

where z/ denotes the exterior normal to G = 'R.^\K. 

Physically speaking, we have that Ki is a sound-soft scatterer, K2 is sound-hard 
and we have an impedance boundary condition on K3. The admissible configuration 
C represents the scattering object for our study, which consists of the impenetrable 
obstacle (i^, 5), the passive medium (cr, g), and the active source/sink {h,H). We note 
the following result that can be proved in a standard way (cf. [27j). 
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Proposition 2.1. Let he an entire solution to the Helmholtz equation Au + k^u = 
in R^. Let us consider an admissible configuration C in R^. Then there exists a unique 
weak solution to the following scattering problem 



\K, 



div(crV^x) + k^qu = —h + div(i7) in 

u = on dKi^ 

aVu ' V — {) on dK2^ 

aVu • V — su — on dK^^ 
^ u — satisfies the radiation condition. 



(2.2) 



If we take u^{x) 



ikx-d 



X G R^, to be the plane wave, then the solution u to 
\2.2\ has exactly the same asymptotic development as that in (1.6), and we denote by 
Uoo{x^d]C)^ X G S^"-*^, the far-field pattern of the scattered field . 

Next, we present a lemma with some key ingredients of the transformation optics, 
the proofs of which are available in [T^l [20l 126] . The main remark is that the definition 
of an admissible configuration is stable under bi-Lipschitz transformations. 

Lemma 2.1. Let and be two Lipschitz domains in R^ and x = F{x) : ft ^ ft be a 
bi-Lipschitz and orientation-preserving mapping. 

Let C = (fi; Ki, ^3, 5, cr, /i, i7) be an admissible configuration in Q. For sim- 
plicity we assume that K and the supports of h and H are contained in ft. We let the 
push- forwarded configuration be defined as 

C^{^-ki,k2.k^,s,a,qXH)^F,{^l-Ki,K2,K:,,s,a,q,h,H) : = 

{F{Sl)-F{Ki), F{K2), F(K3), F,s, F,a, F,q, F,h, F,H). (2.3) 

Notice that Ki is a scatterer, whereas K2 is a regular scatterer and Ks is a Lipschitz- 
regular scatterer. Moreover, 



s{x) = F^s{x) := 



s{x) 



(2.4) 



\det{DrF{x))\ 

where Dt-F denotes the tangential component of the Jacobian matrix of F. We also have 

DF{x) ■ a{x) ■ DF{xY 



a{x) — F^a{x) :- 
q{x) = F^q{x) := 



\det{DF{x))\ 
q{x) 



x=F-^{x) 



(2.5) 



\det{DF{x))\ 

where DF denotes the Jacobian matrix of F. Finally, 

h{x 



x=F-^(x) 



h{x) — F^h{x) :— 
H{x) = F^H{S) :-- 



\det{DF{x))\ 
DF{x) ■ H{x 



\det{DF{x))\ 



x=F-i(J) 
a;=F-l(i) 



(2.6) 



9 



Then u E H^{Vt\K) solves the Helmholtz equation 

V • (aVu) + k'^qu ^ -h^V -H in Vt\K , 

u — Q on dKi, 

aVu • V — {) on 8X2, 

(jVu ' V — su — {) on dK^,, 

if and only if the pull-hack field u — {F~^yu \— uo G H^{Vt\K) solves 



V • {aVu) + k^qu = -h + V • H in n\K, 

5 = on dKi, 

aVu • V — {) on 8X2, 

aVu ' V — ^u — {) on dK^. 

We have made use ofV and V to distinguish the differentiations respectively in x- and 
X- coordinates. 



As a consequence of Lemma |2.1[ one can directly verify that if F : ^ is a 
bi-Lipschitz map with F\q^ — Identity, then the push- forward of an admissible configu- 
ration is again an admissible configuration and 



Uoo{x^d]C) — Uoo{x^d]C)^ for any X, d G ^. (2.7) 



The observation (2.7) is of critical importance for our following general construction 
of the cloaking scheme. Finally, let us point out that the following quantities remain 
unchanged under the push-forward 

I f{Qi)~'\h\^ ^ = 1,2, and [ {a)-^H-H= [jay^H-H. 

With the above preparations, we are ready to present the general construction of our 
proposed cloaks. By a bit abuse of notation, we let from now on fi, D and S be the 

closures of three bounded Lipschitz domains in such that R^\fi is connected and 

00 

H CL)(iri. Let be a compact subset in R^, and be an ^-neighborhood of Kq- 
Both K{) and shall be made precise in the next section. We assume there exists a 
bi-Lipschitz and orientation-preserving mapping such that 

F^ = Identity on R^\f^; F^ on n\K^- F^ on (2.8) 

with 

Fi2)(i^,/2) = E, fP{K,)^D, FPin\K,)^^}\D. (2.9) 

Let, in the physical space, = (Si, XI2, S3, 5, a^, g^, /i, H) be an admissible configuration 
as described in what follows. Set S = XI1UZI2UZI3. The following properties are required. 
First, S C H and h and H are supported in S. Second, we have 

{n\D;a',q^) = {F,)^{n\K,;I,l); a'{x) = I, q^{x) = l for a.e. x G R^\f^, (2.10) 
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and 



(2.11) 



where {K^\K^/2i ^fi Qi) is a lossy layer that is admissible and shall also be made precise 
in the subsequent sections. In the physical space, {^\D; a^, q^) is the cloaking layer, and 
(D\Z;a^^q^) is a layer obtained with specially chosen parameters af and qf. S is the 
cloaked region where the target objects are located, and they include the impenetrable 
obstacle (H, 5), the medium Uf=i Q^) the active source/sink (/i, H). Here, 

we would like to emphasize that in our cloaking scheme, if h and H are zero, namely 
no source/sink is present inside the cloaked region, then (5]\ULl^^;3^^5'^) 
could be arbitrary except for admissibility; and only when there is a source/sink present, 
we would require some generic condition on the cloaked medium Uf=i ^I'l^^ ^ Q^)- We 
finally notice that the configuration inside S do no depend on that is S, 5, h and i7, 
and the coefficients and q^ inside S, do not depend on s. 

The scattering problem corresponding to the cloaking construction described above 
is given by 



div(cr^V'^^) + k'^q'^Ue = -/i + div(i7) 







a^V5e • = 
a^Vu^ • V — su^ = 

Up — satisfies the radiation condition. 



m 



U 



1=1 



on ^Hi, 
on 5112, 
on ^Ha, 



(2.12) 



We assume that the configuration is admissible, therefore there exists a unique solution 
G i7^^^^(M^\ IJ^^^l S/) to (2.12). In the sequel, we let 5^ denote the corresponding 



00)5 



scattering amplitude. 

One of the main results of this work is to show the convergence of (and hence u\ 
and then apply it to derive the approximate cloaking in various settings. To that end, 
in the virtual space, the admissible configuration — {Kf, Kf, ^3, s^, cr^, q^, h^, H^) is 
given by 



that is 



1 = 1,2,3, 



and = KfU K^U C Furthermore, 

a'{x) = I, q%x) = 1 for a.e. x G R^\K„ 



(2.13) 



and 
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We let Ue = {FeYue G HlcO^^\UUi^l) ^^d, by Lemma [2^ it is straightforward 
to verify that is the unique solution to 



Ue = 

a^Vue • = 

(J^VUe • V — S^Ue — 

— satisfies the radiation condition. 



m 



1 = 1 



on OKI 
on J 
on Si^Tf, 



(2.14) 



We notice that (2.12) describes the scattering in the physical space corresponding to 



the cloaking construction, whereas (2.14) describes the scattering in the virtual space. 
Clearly, we have = outside a sufficiently large ball. Hence, in order to study the 
convergence of 5^, it suffices for us to study the convergence of u^, which is the main 
task of Section |3] in the sequel. We note the following peculiar structure of the virtual 
scattering problem (2.14) for our subsequent study. The scattering objects, including 
the passive penetrable medium and impenetrable obstacles Uf=i 5 active 
source/sink are included into the innermost region and then they are enclosed by 

a lossy layer {K^\K^/2i 5 Qi)- When is uniformly small, namely will degenerate 
to a singular point as £: ^ +0, this is the case that has been extensively investigated in 
O El m [20I [211 |23l . As emphasized in the introduction, we shall extend such studies 
to an extremely general setting, especially could be 'partially small', as will be seen in 
our subsequent study. The critical scaling arguments developed in [5l [6l [H |20l [211 [221 [2H1 
cannot be adapted to treating the much more challenging cases of the current study. 



3 Analysis in the virtual space 

This section is devoted to the convergence analysis of the virtual scattering problem 
(2.14). However, our study shall be conducted in a much more general setting than is 
needed for the cloaking purpose. 

Let Kq be a scatterer in R^. Let us denote by d : [0, +oc) the distance 

function from Kq defined as follows 



d{x) = dist(x, Kq) for any x E 



We assume that there exists a Lipschitz function d : R^ [0, +oc) such that the 
following properties are satisfied. 

First, there exist constants a and 6, < a < 1 < 6, such that 



ad{x) < d{x) < hd{x) for any x G 



For any s > 0, let us call = {x E R : d{x) < e}. For some constants ^0 > 0, 
p > 2, Ci > and i? > 0, we require that for any 5, < 5 < ^o, i^^ C Br^ W^\K^ is 
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connected and 



Alv{Br+^\k,) < Ci\\u\\j^i^Bn+i\K,) for any u G H^{Br^i\K^). 



We notice that, under these assumptions, K^, for any < 5 < Sq, is a regular 
scatterer. A simple sufficient condition for these assumptions to hold is that Kq is a 
compact convex set. In fact, clearly we have that M.^\Ko is connected. Then, we can 
take d — d OY the distance from with respect to any norm on M^, not only with 
respect to the Euclidean one. Then, for any 6: > 0, = B^(Ko)^ clearly with respect 
to the chosen norm. For any 6 > 0, we have that B^{K) is still a convex set, therefore 
MJ^\B^{Kq) is connected and satisfies a cone condition with a cone independent on e. 
Hence also the other required properties are satisfied, for instance using [1, Theorem 5.4]. 

Another sufficient condition is that is a Lipschitz scatterer, see [30, Section 4]. 
Finally, may be the union of a finite number of pairwise disjoint compact convex sets 
and Lipschitz scatterers. 

For any fixed Q < e < e^^ let us assume that {A^f , A^f' ^ 
admissible configuration. 

Moreover, we require that one of the following two assumptions holds. 

Assumption 1. {Kf , Kf, ^f, 5^, cr^, g^, /i^, is an admissible configuration with the 
following properties. First, = U U C Second, the behavior of and 

is different in the following three regions, K^j^K^ ^ K^\K^i2 and 'R.^\K^. 

a) We have that a^{x) = / and q^{x) = 1 for almost every x G MJ^\K^. 

b) There exist a continuous nondecreasing function uji : (0,6o] (0, +oc), such that 
lim5^o+ ^il"^) — positive constants £"1, A, and E2 such that for almost 

any x G Ke\K^/2 



< ql{x) < Ei{cji{£))-^ and < E[{cji{e))-^ < ql{x). 



(3.1) 



Furthermore 



(y'{x)i • ^ < A||^f for any ^ e and for a.e. x £ Ke\K^/2 



(3.2) 



and 



^cr^(x)vJ(x) • Vd(x) < E2 for a.e. x G K^\K^/2. 



(3.3) 



Finally we require that 




(3.4) 



c) There exists a positive constant E3 such that 



< sup qf < E3 inf g|. 
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We assume that = for almost every x E \K^/2- Moreover, we have 

that 

{k^q'^y^lh'l^ <E3 and / {a'Y^H' - W < E^. (3.5) 

Here the constants Ei, E[^ A, E2 and E^ do not depend on e. 

We remark that the previously stated assumption, essentially because q2 > in 



K^\K^, implies that the corresponding problem (2.1) has only the trivial solution with 



out any further assumption on the obstacles Kf, Kf? ^3 the coefficients and a^. 



Assumption 2. Assumption [2] coincides with Assumption [l] except for point c) which is 
replaced by the following. 

c) We assume that = and that h^{x) = for almost every x G R^\K^/2' Moreover, 
there exists a positive constant Es such that 

/ {k'ql)-'m'<Es. (3.6) 

Notice that the above condition means that q2 > almost everywhere on the set 
{h^ ^ 0} and that the integral is actually performed on such a set. 

Remark 3.1. Let us point out here the main differences between these two assumptions. 
In Assumption [2] we do not impose any condition on the coefficients and scatterers inside 



K^/2 except for the indirect condition in (3.6). In particular, if there are no sources, then 



no condition at all is required for the configuration inside K^/2 except for its admissibility. 
The drawback is that the kind of sources we can allow in Assumption [2] is more limited 
than those considered in Assumption [TJ 

Let us finally notice that, only for Assumption [l] we may replace ( |3.3| ) with the 
following weaker condition 

lim \ [ a'Vd-Vd^O. (3.7) 

We shall keep fixed throughout this section two functions (/>!, (f)2 ^ C^(R) such that 
they are nondecreasing and is identically equal to on a neighborhood of (— oc, 3/4] 
and identically equal to 1 on a neighborhood of [1, +oc), whereas 02 is identically equal 
to on a neighborhood of (— oc,l/2] and identically equal to 1 on a neighborhood 
of [3/4, +00). Also for any symmetric positive definite matrix a, we define ^/o^ the 
symmetric positive definite matrix such that ^/a • ^/a = a. 

Let us also fix < 5^?, < ^o? ^ ^ a decreasing sequence of positive numbers 
such that lim^ £n — 0. For simplicity we denote 

Kn i^6n and Hn K^^/2 and ^36^/4^ 
We begin with the following result where a Mosco-type convergence is proved. 
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Lemma 3.1. Under the previous assumptions^ the following two properties hold. 

i) Let us consider a subsequence 8^ and a sequence ui G H^{BiiJ^i\Hn^), I G N. Let 
us consider (ui, \/~a^Vui) G L'^{Brj^i^M.^^^) by extending ui and \/~o^Vui to in 
Hni- If {ui, ^/a^Vui{l — XJni)) converges weakly to (u^V) in L'^^Br^i^M.^^'^) as 
I OC; then u G II^{Br^i\Ko) and V = Vu in Br^i\Ko. 

a) For any (f G II^{Bji^i\Kq) such that (f is bounded, there exists a sequence (fn ^ 
II^{B Rj^i), n G such that Lpn — ^in Jn and {(fn, Vo^^y^n) converges strongly to 
{(f,\/ip) in L^(5i^+i, R^^-'^) as n ^ oo, where (f and Vip are extended to in Kq. 

Proof We begin by proving[i]). Let D be any open subset of Br^i such that DDKo = 9. 
There exists I G N such that for any / > I we have that, in {ui, \/a^i\/ui{l — XJni)) ~ 
{ui,Vui). Therefore V = Vu in D and we may easily conclude that u G II^{Bji^i\Ko) 
and V = \/u in Bji-^i\Kq. Moreover, v and V are almost everywhere in Kq. 

For what concerns , we define Xn — and (fn — Xn^- We immediately 

obtain that is identically in and that c/:?^ = in BR^i\Kn. Clearly (fn ^ 
II^{B Rj^i) and, as n ^ cx), converges to in L'^{Bji^i). Then we need to estimate 



The first term of the right-hand side clearly goes to as n ^ oc. About the second 
term, we may notice that 

WV^Vifnf < 2 (||(V^V(/^)Xnf + ||(V^VXn)^f 

and that 

^ ///7/ 

VXn Md/en) — . 

The conclusion immediately follows by ( |3.2| ) and ( |3.7[ ). □ 
Let u'^ be an entire solution to the Helmholtz equation Au + k^u = in R^. We 



consider the scattering problem (2.2) with the configuration {Ki, K^, a, H} 
replaced by {A^f , 5 s^^a^^q^^h^^ H^} and we denote by Un its solution. 
We have the following theorem. 

Theorem 3.2. Under the previous assumptions, Vn = Un{l — XHn) converges to a func- 
tion u strongly in L'^{Br) for any r > 0, with u solving 

Au + k^u^Q in R^\Ko, 

u^u'^u' in R^\Ko, .^.x 



lim,^+oo r(^-i)/2 - iku') = 



r — \\x\ 



and u — {) in Kq. 
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Proof. We develop the proof when Assumption [T] holds true. At suitable points of the 
proof, we shall describe the needed modification if Assumption [2] is used instead. 
Let us begin by assuming the following further condition 



Ml^Br^^) < C for any n G N. 



(3.9) 



Let us assume that, up to a subsequence, we have that Vn converges weakly to u in 
L?{B Rj^i) as n ^ oc. We immediately infer that u — in K^. 

We notice that Un satisfies the Helmholtz equation /\u + k^u = in W^\Br. By 
standard regularity estimates, we may infer that there exists a constant Ci, depending 
on the constant C in (|3.9[) and i?, such that for i?i = i? + (1/2) we have 



Wn\\c^{dBR^) + W^UnWc^^dBR^) < Ci for any n G N. 



By the weak formulation of (2.2), we infer that 

Jbr^\k^ 



(3.10) 



where we recall that K"" = Kf U U 



Clearly the exists a constant C2, depending on Ci and R only, such that 



Jobr^ 



< C2 for any n G N. 



We need that {{vm Vo^'^Vni^ — XJn))}neN is uniformly bounded in L^(Si^^ , R^+-^), 
that is 

\\{vn, \^Vvn{l - XJn)) II ^ ^^r any n G N. (3.11) 

Notice that we are assuming as usual that y/a^Vvn is extended to zero in Hn. If this is 
the case, up to a subsequence, (i;^, \/^\/Vn(l — XJn)) converges weakly in L'^{Br^ , R^+i) 
to (v^V) as n ^ oc. Clearly v — and, by property of Lemma [3^1} we have that 
u G H^{Br^\Ko) and V = Vu in Br^\Ko. Moreover, u and V are zero in Kq. 

Let us prove such a uniform boundedness. Here the proof is slightly different de- 
pending on the Assumption used. We begin with Assumption [T] 

Obviously 



0< / k'^q^'Kl'^ 

lBR^\Kn 



ifi I ^ k C . 



BR^\Kn 
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We split into the real and imaginary parts and, by the properties of and 5^, we obtain 
that, inside Kn^ lHn\K^ ^2 I'^riP and ixn\Hn ^2 I'^riP are not negative. Furthermore, 



0< / k\^\un\^+ [ eq^\un 



for any positive constant c, which we shall fix later. Hence, we have that both integrals 
Ih \K'^ ^^^2 I'^nP and A;^g2 I'^nP are bounded by the following quantity 

An = Es + — +C [ (J^'VUn • VU^ + 2C2. 

We conclude that 

I \un\^ <{An/E[)uji{en) and f \un\^<An( inf q^] . (3.12) 



Hence 



/ k^'q^lunl^ <{Ei/E[)An and / k\^\unf < E^A^. 



We deduce that 



1, 

k^C^ + {Ei/E[)An + EsAn + ^ + + ^"V^, • + ^2 

Z Z ZCi Z JHn\K^ 



for any positive constant ci, that we shall fix later. Setting E = {Ei/E[) +Es + l/2 and 
C = /c^C^ + £^3/2 + (72, we choose c and ci such that Ec + ci/2 = 1/2 and we obtain 
that 

/ a'^Vun • Vi^ < 2(7 + 2^^i + ^Js/ci 

where £'1 = £'3 + £'3/c + 2C2. 

We conclude that there exists a constant Ci such that 



thus in particular 



/ a^Vun • V?iri < Ci for any n G N 



/ cr^'Vi^n • Vi^ < (5i foranynGN (3.13) 

^ Bji-^\Jn 
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and the required uniform boundedness property is achieved. 

Moreover, we can also infer that there exists a constant A such that An < A for any 
n G N and that 



/ knP < iA/E[)uJi{en) for any n G N. (3.14) 

J Kn\Hn 



Finahy, there exists a constant C2 such that 



and in particular such that 



/ k^\q''\\un\^<C2 for any n G N 



/ < ^2 for any n G N. 

^ Bji^ \Hn 



(3.15) 



For what concerns Assumption [2j the reasoning may be changed in the following 
way. We have that 

0< / k\^\un\^+f eq^\un\^ K-^ + l f k\^\un\^ + C2. 

Therefore, we obtain that 



< 



/ k^q^W 

J Kn\Hn 



We easily conclude that (3.14) and (3.15) hold true also in this case. 



Let us now define, for any n G N, Xn = 4^2{d/£n)- We apply a Caccioppoli-type 
inequality. Namely, we have that 

JBr^XHu JBr^XHu JdBR^ 



Therefore, 



Br^ \Hn 



< 



C2 + C2 + 



2 JBR^XHn 



Br-^ \Hn 



that is 



X^^'VUn ■ < 2C2 + 2C2 + ^ / \Un\\4>2{d/en)Ycj''yd ■ Vd. 

Br^\H„ Jk„\H„ 



By (3.3), we conclude that (3.13) also holds true. 
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Let us fix (/p G H^{Br^i\Ko) whose support is contained in Br-^ and such that (f 



is bounded. By property p| of Lemma 3.1, we can construct for any n G N a function 



(fn ^ H^{Brj^i)^ whose support is contained in Br^ and such that = in J^, with 
{(fni Vo^^^n) converging strongly to in L^(Si^+i, M^+-^) as n ^ oc, where (p 

and \/(f are extended to in i^o- 
We have that 

/ a'^VVn • V(Pn - / k^q'^Vn^n = 0. 

^Br^ Jbr^ 

It is easy to show that 



hm 



Jbr^ Jbr^i\Ko 



We observe that 



/ k^q'^VnPn = / k^VnPn + / k^q'^UnPn = 

JBr^ J Bpi^l\Kn J Kn\Hn 

J Br^I J Kn\Hn 



Obviously 

/ ^ri^^BR^,\Kn = / 

^Br^I J Br^i\Ko 

Finally, since (fn is uniformly bounded and k'^\q'^\\un\'^ is uniformly bounded, and 



by (|3.4|), we are able to pass to the limit and prove that 

Vu • V(f — k^u(f — 0. 



^i?+l\^0 

By density, we have this last equation holds for any (f G H^{Bji^i\Kq) whose support 
is compactly contained in Br-^ . 

By Lemma 3.1 in [37j, we have that, up to a subsequence, Un converges to a function 
u uniformly on compact subsets of M.^\Br, with u solving 

Au + k^u = in R^\B^ 

u = u' + u' in R^\B^ 

limr^+^A^-^y^{^-iku') =0 r= ||x||. 



Then we immediately conclude that u is the unique solution to (3.8). We now 
show that the whole sequence Vn converges to u strongly in L?{Br) for any r > 0. By 



uniqueness of the solution to (3.8), we have that the whole sequence Vn converges to 
u weakly in L'^{Bji^i/2) and uniformly on compact subsets of M.^\Br. Moreover, by 
standard regularity estimates, we may also assume that Vn converges to u strongly in 
on any compact subset of Bji^i\Ko. 
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We recall that there exist constants p > 2 and C4 > such that for any n G N we 
have 

\H\LP(Bn^\Kn) < C4\\u\\Hi(Bn^\Kn) ^ ^ H^{BR^\Kn). (3.16) 

Let us fix a positive constant /iq, < /iq < minj^o, 1/4}. For any < h < ho, let 
Ah = Kh U {BR-^\Bji^_h) and S/^ = Bn-^Ah. We notice that B^ is compactly contained 
in Bji^\Ko and that lim/j^o+ |^/i\^o| — 0. 

Therefore, for any /i, < /i < /iq, there exists n G N such that for any n >n we have 
En ^ h and 

ll^n - ^/||l2(Bh^) < \\u\\L^Ah\Ko) + ll^n||L2(K^\i^^) + \\Un\\ L^{Ah\Kn) + ll'^ri " '^IIl2(B^). 



By (3.16) and (3.11), we infer that 

\lHBr^) < \ML^Ah\Ko)HUn\\L^Kn\Hn)+C3C4\Ah\Kn\^^~^^^^ 



By ( |3.14[ ), we conclude that for any n > n 

\(p-2)/(2p) 



U\ 



L^A^) < \H\lha^\Ko) + J {A/{k^E[))uJi{Sn) + CsC4\Ah\Ko 



We fix 5 > 0. There exist /i, < /i < /iq, and n E N such that for any n > n we have 
Sji ^ h and 

There exists h>n such that for any n>h we have 

\\Un - '^||l2(B^) < 5/2, 

therefore we have proved that Vn converges to u strongly in LP'{Br-^) and consequently 
strongly in L'^(Br) for any r > 0. 

In order to conclude the proof, we have to show that (3.9) holds true. Let = 
lbn||L2(Bi^+i)- By contradiction, let us assume that lim^^oo = +oc, possibly by 
passing to a subsequence. Let us consider Wn = Un/dn- We have that 

-Xifn)llL2(BH+l) = 1- 

Therefore it;^(l — XHn) converges to a function w strongly in L? on any compact 
subset of R^. The function w satisfies 

Alt; + fc^it; = in R^\Ko 

Clearly we also have that — 1- 

We have that ||'^n/^n|li,2(^^^^\;g^), n G N, is uniformly bounded. Therefore, again 
up to a subsequence, uf^/an converges, as n ^ oc, to a function w strongly in on any 
compact subset of R^\Br. Such a function w satisfies 



Aw -^k^w^O in \Br 

lim,^+oor(^-')/2(^-zfcri;) =0 r=||x||. 
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Since Wn = u'^ /an -\- u^/an, we may immediately conclude that, outside Br, we have 
w — w. That is w solves 

Aw + k'^w = in R^\Ko 

Vit;-zy = on^Ko (3.19) 

Imir^+ooT^^-^^'^ - ikw) =0 r= ||x||. 

By uniqueness, and since R^\i^o is connected, we may conclude that w is identically 
zero, which leads to a contradiction since ||'?^||L2(Bi^^i) should be equal to 1. □ 

We summarize the results of this section in the following theorem. 

Theorem 3.3. Under the previous assumptions, for any e, < e < e^), let us fix 
af = erf (^); X G K^\K^i2, an N X N symmetric matrix whose entries are real-valued 
measurable functions such that, for some X^, < < 1, we have 

^eUf <^'{x)^'^ for any ^eR^ and for a. e. x e K,\K,/2. 

Let qf = + iq2 = Qi(x), x G K^\K^/2^ ^ complex- valued bounded measurable 
function, with real and imaginary parts qf and ^ respectively, such that, for some X^, 
^ < X^ <\, we have 

^i(^) ^ X^ and ^l(x) > for a.e. x G K^\K^i2. 

We further require that af and qf satisfy Assumption ^p), for a given function uji 
and given constants Ei, E[, A, and E2. 

Then let {Kf,K2^K^,s^,a^,q^,h^,H^} be any admissible configuration such that 

= a I and q^ = qf in K^\K^/2 

and satisfies, for a given constant Es, either Assumption [T] or Assumption [2} Notice 
that if Assumption^ is used, we may replace (3.3) with ( |3.7| ). 

For any d G let u'{x) = e^^^'^, x G W^. Let Ue be the solution to ([2^ with 

the configuration {A^i, A^2, ^3, <§, a, /i, H}. replaced by {A^f , A^f? ^ H^}^ 

and u be the solution to (3.8). Let and u^o be the scattering amplitudes of the 
corresponding scattered waves, respectively. 

Then there exists a function uj : (0,£o] — ^ (0, +oc] such that lim5^Q+ 6^(5) = 0^ 
depending on Kq, d, R, af , qf, and E^ only, and a constant C depending on R only, 
such that for any e, < e < eo, we have 

and 

\\ulo - '^ooIIl-(s^-i) < Cu{e). (3.21) 
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4 Full and partial cloaks in the physical space 



In this section, we shall apply Theorem 3.3 to the constructions of full and partial 



invisibility cloaks in the physical space. We shall stick to the terminologies intro- 
duced in Section [2] for the general construction of the cloaks. Specifically, we shall let 
{Kf , K|, Kf, 5^, cr^, g^, /i^, H^} and {Hi, 112, II3, 5, a^, g^, /i, H} be the scattering objects 
respectively in the virtual and physical spaces. Let be the blow-up transformation 
from the virtual space to the physical space. 

4.1 Regularized full invisibility cloak 

Suppose D is the closure of a bounded convex domain in that contains the origin. 

o 

Let ft be the closure of a bounded Lipschitz domain in R such that D Cfl and ^\D is 
connected. 

For any s > 0, let = {ex;x G D}. It is supposed that for any < e < 1 there 
exists a (uniformly) bi-Lipschitz and orientation-preserving map, 

fP : n\De n\D, FPldQ = Identity. (4.1) 
That is, fJ^^ blows up to D within ft. Let 

fP{x) = -, xeD„ (4.2) 



e 



and 



(Identity on R^\n, 
fP on n\D,, (4.3) 



1 



fP on D,. 



Formally, we set 
Furthermore, we set 

and 
with 



E = I?i/2, Ke = D„ K,/2 = D,/2. (4.4) 

al^{FP)J, g^ = (FP),l mn\D, (4.5) 

iD\D,/,;af,qf) = {FMKe\K,,^;at .qf), (4.6) 

Cjf = Cl(x)£2, gf = (C2(x) + iC3(x))£-^+\ (4.7) 



In (4.7), ci(a;) is a symmetric-matrix valued measurable function, and C2{x),c^{x) are 



bounded real valued measurable function such that 
Aolief < ci(x)e-C < Aolief , Ao < C2(x),C3(x) < Ao for a.e. x G i^A^./2, (4i 
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where Aq and Aq are two positive constants independent of e. {D\Di/2]d^i ^qf) is the 
chosen lossy layer for our cloaking scheme. By Lemma [2. If it is straightforward to verify 
that 

af (x) = ci{£x)s^, qj{x) = {c2{£x) + ics{£x))£, X E D\Di/2. (4.9) 

We would like to emphasize that ( |4.9[ ) is only a particular choice for illustration of our 
general results in Section [sj and there are more choices as long as (K^\K^/2i '^f ^ Qi) is 
such chosen that Assumption 1, b) is satisfied. 

Next, we first consider the cloaking of passive objects by assuming —h + div(i7) = 
in the physical space. Let (Si, XI2, S3, s) be an admissible obstacle located inside D1/2, 
and {Di/2\[J^=i^i'-,^a^Qa) be an arbitrary regular medium. Let (K^^K'^^K^^s^) and 
{Ks/2i^a^Qa) be the corresponding virtual images in the virtual space. Also, we set 

= Uf=i ^f- Hence, in the physical space 



1. 1 



m 



pN 



%ql in n\D, 

in D\ 
Ga.qa in 



af,qf in D\Di/2, 



(4.10) 
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Proposition 4.1. Let {Si, XI2, S3, 5, a^, g^} and {K^,K^,K^,s^,a^,q^} be as described 
above, andu^ and be the corresponding scattering waves [cf. ( 2.12[ ) and (2.14)). We 
also let 5^ and u^^ be the scattering amplitudes of, respectively, and u^. 

Then there exist > and a function uj : (0, Sq] ^ (0, +oc] with lim5^o+ ^{^) — 0; 
which is independent ofaaiQa ci'^d Tii,T^2,^3, s, and of the impinging direction d, such 
that for any £ < £{) 

= < ^(^)- (4.11) 



Proof. Let O be the origin. For any x E M^, let x E dD be the point lying on the line 
passing through O and x. Set 

d{x) = dist(x, 0)~^d{x) for any x E R^. 

By taking Kq := {O}, it is easily seen that = and = Since Kq has zero 

: 

□ 



capacity, we know that the solution to ( |3.8[ ) is given hy u — 
and 



0. One readily has (4.11) by Theorem 3.3 



Consequently, 



Proposition |4.1| indicates that one would have an approximate full invisibility cloak 



for the construction (4.1 )-(4.10). The essential point in Proposition 4.1 is that a single 



point has zero capacity. We recall that H^{D) — H^{D\Kq) for any open set D and any 
compact Kq C D with zero capacity. We would like to emphasize that by following the 



same spirit, and using Theorem |3.3| one could have more approximate full invisibility 
cloaks. For example, in R^, a line segment is also of zero capacity, and hence one 
could achieve an approximate full cloak by blowing up a 'line-segment-like' region in R^, 
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namely Kq is a line segment; or by blowing up a finite collection of 'point-like' and 'line- 
segment-like' regions. In Section [5| we shall give more discussion on how to construct 
an approximate full cloak by blowing up a 'line-segment-like' region. 
Finally, we consider the cloaking of active contents. 



Proposition 4.2. Under the same assumptions of Proposition |4.1t we further assume 
that the source/sink term satisfies the following. We let H — {) and h G L^(R^\5]) be 
such that supp(h) C Di/2\^- Moreover, we require that 

^Qa ^ Aq > on supp(h)^ (4.12) 

where Xq is a constant. 

Then there exist £q > and a function uj : (0,£:o] (0, +oc] with lim^^o+^l"^) — 
0; which is independent of aa^^qa^^^a\{h=o} 5]i,5]2,S3^ s, and of the impinging 
direction d, such that for any s < 

= < ^(^)- (4.13) 

Proof. It is straightforward to verify that 



/2 



(4.14) 

Hence, Assumption 2, c) is satisfied. The proposition can be proved following a com- 
pletely similar argument as that for Proposition |4. 1 1 by using Theorem |3.3[ □ 



Remark 4.1. Proposition |4.2| indicates that in addition to passive mediums, the con- 



struction (4.1)-(4.10) is also capable of nearly cloaking active contents. We need that 
at the place where the active source is located, the medium must be absorbing. It is 
recalled that the Helmholtz equation can also be used to describe the electromagnetic 
phenomena. In such a case, h denotes an electric current density, whereas S^^a denotes 
the conductivity. At the place where '<sqa = 0, one must have /i = since there the 



medium is non-conducting. Hence, (4.14) is a reasonable physical condition. Moreover, 
we would like to emphasize that one can nearly cloak a more general source term of the 
form — /i + div(i7) by using Assumption 1, c) by imposing certain physical conditions on 
a a and 5a- 

4.2 Regularized partial invisibility cloak 

Our construction of partial cloaking devices will rely on blowing up 'partially' small 
regions in the virtual space. We first present our study in the virtual space. Let 

Ko := {-a < xi < a} X {x2 = 0} in (4.15) 

and 

Ko := {-a < XI < a} X {-b < X2 < b} x {xs = 0} in R^. (4.16) 
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We note that v = (0, 1) in 2D and v = (0, 0, 1) in 3D for K^. Let < r < 1 and 
define 

AA^ {0 G S^-^ : \v'e\<T}. (4.17) 



Next, we consider the scattering problem (3.8) with given in (4.15) and (4.16), 



which is known as the screen problem in the hterature. 



Proposition 4.3. Let Kq be given in (|4.15|) and (|4.16|)^ and Mr be given in (4.17). 

ikx-d^ T/ien there exists a 



Kq) be the solution to (3.8) with u^{x) — e 



Let u e Hl^ 

constant C , depending only on a, h and k, such that 

\uoq[x^ d)\ < Cr for x G S^"-*- and d G A/^, 

and 

\u^{x^ d)\ < Cr for x G A/V and d G Sat-i- 
Proof. By the weh-posedness of the forward scattering problem, we know 



, du 



«oo(-,rf)||Loo(siV-i) < C||^||i2(^„) = Ck\v ■ d\{H''-\Ko)) 



.1/2 



(4.18) 
(4.19) 

(4.20) 



from which one readily has (4.18) by recalhng the definition of Mr ■ (4.19) can be obtained 

(4.21) 



from (4.18) by using the following reciprocity relation (see, e.g., [9J) 

Uoo{—d, —x) — Uoo{x, d) for x G and d G S^"-*^. 



The proof is completed. 



□ 



Now, the construction of a partial cloak shall be based on the use of Theorem |3 . 3| and 
Proposition |4.3| similar to the one for the full cloaks in Section [4T] by following the next 
three steps. First, one chooses K^, an ^-neighborhood of Kq, and a blow-up transforma- 
tion F^, and through the push- forwards, one constructs the cloaking layer {^\D; a^, q^). 
Second, according to Assumption 1, b), a compatible lossy layer [K^\K^i2]cr^ ^q^) is 
chosen in the virtual space, and then by the push-forwards, one would have the cor- 
responding lossy layer in the physical space. Finally, according to Assumption 1 or 2, 
c), one can determine the admissible media, obstacles, or sources that can be partially 
cloaked. Next, let us consider a simple 2D example by letting = [—1—6, l-\-s\ x [— s], 
which is an ^-neighborhood oi — [—1^1] x {0}. If a uniformly blow-up as the one in 
(4.2) is used, then will be transformed into [— 1/s — 1, l/s+ 1] x [— 1, 1]. In such away, 
one would have a partial cloaking device of a very large size depending on . For prac- 
tical considerations, in the next section, we shall develop an assembled-by-components 
technique on blowing up a partially small region in constructing a regularized partial 
cloak of compact size. 
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5 ABC geometry and concatenated construction of three 
specific cloaks 

In this section, we shall first discuss the assembled-by-components (ABC) geometry in 
the virtual and physical spaces, and then give the blow-up construction of three specific 
cloaks based on concatenating individual components. We would like to emphasize that 
we would not appeal for a most general study in this aspect, and instead we shall present 
our study based on these three specific examples. However, the technique developed 
could be straightforwardly extended to devising many other cloaks. 

5.1 ABC geometry in the virtual and physical spaces 

Define for x = (xi, . . . ^xn) G M^, 

\x\w,oo '= m3x{\wiXi\}iL^, (5.1) 

and for 1 < p < +oo, 

\^\w,p •= I X] l^ml^ I , (5.2) 

where w = {wi^ . . . , wn) G with < wi < 1 denotes a weight. For any r > we 
define the (i(;,P)-ball and sphere, respectively, as follows 

n^p)^ := {x G : \x\^^p < r} and S^P], := dn^P^, 1 < p < oo. (5.3) 

Two semi-(tf;, P)-balls can be defined by 

f^Sj ■■= e ^(^1. : Xk>0} and f^^^^- := {x G ^^^^ : < 0}. 

(v) 

^w!r with p = l,2,oc will be the three base geometries for our subsequent construc- 
tions. Next, we show that by concatenating the base geometries, one could obtain more 
practical geometries in both the virtual and physical spaces. 

In the sequel, we denote x^ := {xi^ . . . ^xj^-i^x^^i^ . . . ^xjsf) G M^"-"^ where the su- 
perscript k indicates that the fc-th component x^ is dropped. Thus N — 1 dimensional 
{w^^F)-h3lls and semi-balls can be defined, respectively, by 

n^Pl := {\x\, „ < r} and f^^f ^. := {x^ G O^^^ : xj ^ 0}, 

where j is different from k. 

5.1.1 2D ABC example C 

We start with a simple 2D example. Let W{) := [1,1] and 

Cl = ^^ilvTi + [-«> 0] , Cr = + [a, 0] (5.4) 
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and 

Cm {—a < xi < a} X {—r < X2 < r}, (5.5) 

where r and a are positive constants, that we shah specify in the fohowing. By assembhng 
the three components, we obtain 

C[^o,r,a] = Cl U Cm U Cr. (5.6) 




(a) (b) (c) 

Figure 1: Two dimensional ABC geometry C depicted in both the virtual space (between 
innermost and outermost boundaries) and the physical space (between intermediate and 
outermost boundaries). From left to right: and cloaks, respectively. 

If one takes p — p' — 2^ then C[wQ^r,a] is a 2D capsule-^ if one takes p — oo and 
y = 1, then C[wQ^r,a] gives a 2D nail] and \i p — p^ — oc, then C[wQ^r,a] is a rectangle] 
see Figure [l] for these three constructions. In our subsequent construction of the partial 
cloaking device, we would take C[wQ^e,a] the 'partially' small region in the virtual 
space, whereas C[^q ,p2,a]\C[^Q,ri,a] is the cloaking region in the physical space. Here 
r2 > ri > e with r/ ^ 1, / = 1,2 and a is a free parameter ranges from £ to 1 that 
determines the apertures of the cloaking device. It can be easily seen that if a = 1, then 
C[t(;o,£,i] is 'partially' small since it is only small in one dimension and of regular size in 
another dimension, whereas ii a — e then it is uniformly small since, along both the xi- 
and X2-dimensions, the region is small. 

5.1.2 3D ABC example D 

We next consider a 3D example. Let W{) \— [1,1,1] and define three detached components 

Dl = fiilvTi + [-a, 0, 0] , = f2i^j;+ + [a, 0, 0] (5.7) 

and 

Dm = {-a < xi < a} X {x^ e (5.8) 

where r and a are positive constants, that we shall specify in the following. By assembling 
the three components as shown in Figure [2j we obtain the slender cloak 

^[w,r,a] = Dl U Dm U T>R. (5.9) 
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(a) Detached components in Type D cloak. 




(b) Slender cloak with assembled components. 
Figure 2: Schematic illustration for the construction of Type D cloaks. 
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For more concrete examples of this slender cloak construction, let wq = [1,1,1], 
a = 1. We take r = e and r = 1, respectively, to represent the geometries in the virtual 
and physical spaces. If p = ^ = 2, then D[^q ij is a 3D capsule; if p = oc and = 1, 
then D[^Q 1] is a 3D nail; and if p = ^ = oc, then Dj^^ ,^ is a rectangular-prism; see 
Figure [31 for illustration. In our subsequent discussion, Dj^^^^ will be the region in the 
virtual space whereas Dj^^ ^ will be the region in the physical space. As can be easily 
seen, in the virtual space, Dj^^ ^^^j is small along the X2-, xa-dimensions while it is large 
along the xi-dimension. For our subsequent study, we are also interested in the case 
when the region in the virtual space is only small in the X3-dimension while it is large 
in the xi-, X2-dimensions. It is directly verified that D[^^ i^^j with = [1,5, 1] satisfies 
such requirement. However, it can also be verified that for such w^, in the physical space 
D[^^ 1 1] is very large along the X2-dimension, which is actually of size In order 

to construct a more practical partial cloaking device with its size independent of the 
regularization parameter 6:, we would like to further develop the ABC geometry in the 
following. It is emphasized again that this is our major motivation in developing the 
ABC geometry. 




Figure 3: Three dimensional ABC geometry D depicted in both the virtual space 
(between innermost and outermost boundaries) and the physical space (between inter- 
mediate and outermost boundaries). From left to right: and l^ cloaks, respectively. 



5.1.3 3D ABC example E 

We next consider a specific example for further developing the ABC geometry (see 
Figure [i]). Let 

Eo = {-a < xi < a} X {-b <X2<b}x {-r < xs < r}, (5.10) 

where r = s or r ^ 1, and a, 6 are positive constants. For x = [xi,X2,xs] G M^, let 
w := [1, 1, 1] and 

E^ ={x : x^ E ± [a, 0,0], -6 < X2 < 6}, 

(5.11) 

Ef={x: Gf^^t; ±[0,6,0],-a<xi <a}. 
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(a) Detached components in Type E cloak. 




(b) Flat cloak with assembled components. 
Figure 4: Schematic illustration for the construction of Type E cloaks. 
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and 



={x + [a, ±6, 0] : X G ^^^^ with xi > and ± X2 > 0}, 

, . (5.12) 
E4 ={x + [-a, ±6, 0] : X G with xi < and ± X2 > 0}. 



Let 



^[wo,r,a,b] 



Eou|^|jEfj. (5.13) 



E is a rescue cushion composed of nine components: the centered component Eg which 
is a cuboid of dimensions 2a, 2b and 2r; four components E^2 attached to the side faces 
and four components E^4 assembled to the four corners. In our subsequent study, we 
would take r = 6 as the geometry in the virtual space while r ^ 1 as the geometry in 
the physical space. 

5.2 Blow-up construction and concatenation 

We shall be concerned with two basic types of blow-up transformations. Henceforth, let 
r2 > ri > e > 0. Set 

F,(x) = (A + B\x\^^p) ^— , (5.14) 

\X\w,p 

where 

A=^r„ 3=^:^. (5.15) 

r2 - £ r2 - £ 

It is verified directly that blows up ^w^e to ^w\i within ^w^r2 5 namely 

Fe (n^Pl^Xn^^l^ = = Identity. (5.16) 



For A and B given in (5.15), we let 



\ XI, 1<1<N-1 ^ ^ 

and it can be verified that 

Ge {QrAQe) = QrAQn, G.|n., = Identity, (5.18) 

where 

Qr :^{\xi\ <ai: 1 < I < N - 1} x {\xn\ < r}, 
:={|x/| <ai: 1 < / < TV - 1} x {|x7v| = r}. 



Now, by using blow-up transformations of type ( |5.14 ) or type ( |5.17 ), and concate 



nating each sub-component, it is straightforward to see that there exists a bi-Lipschitz 
and orientation-preserving map ^i, such that for the region constructed in (|5.6[) 



^1 (C[^o,r2,a]\C[^o,e,a]) = C[^Q^^2,a]\C[^o,ri,a], ^il^Cf ] = Identity. (5.19) 
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In a similar manner, one can show that for D and E, respectively, in (5.9) and (5.13), 
there exist and ^ such that 



^1 (D[^,^2,a]\D[^,e,a]) = D[^^^2,a]\D[^,ri,a], ^1 laDf^,,^,^] ^ ^entity, (5.20) 



and 



^ (E[^o,r2,a,6]\E[^o,5,a,6]) = '^[wo,r2,a,b]\^[wo,ri,aM^ ^l\dE.^^^^r2,a,b] = ^entity. (5.21) 

5.3 Partial cloaking devices 



With the preparations in Sections |5 . 1 1 and [5^ we shall present our cloaking construction 
in this section. We start with a simple 2D case. Let be given in (5.19) and set 



'identity for x G R^^Cf^Q^^^,^]^ 

^l(x) for X G C[^Q^^2,a]\C[^o,6,a], 

^^s{x) for X G C[^,,e,a] = Cl U Cm U Ci^, 



(5.22) 



where 



xi 

X2 



ri 

£ 



{xi + a) — a 



ri 



X2 



^ < 



Xi 

ri 
e 
ri 

£ 



X2 J 

(xi — a) + a 



ri 



X2 



for X G Cl, 
for X G Cm, 
for X G Ci^. 



(5.23) 



Clearly, we have that ^ : C[^q ,^2,a] ^ C 
and ^|aC[^Q ^2 a] ~ Identity. Next, let 



^1 is bi-Lipschitz and orientation-preserving 



ci{x) 

C2(x)5 



2 I ^ (^) = (C3(x)+ic4(x))5 



where Q, 1 < / < 4 are positive functions satisfying 

Ao < Ci{x) < Ao for a.e. X G C^wo^eM\^[wo.e/2,a]^ 

with Ao and Aq two positive constants independent of e. Then we let 
We further let 



(5.24) 



(5.25) 



(5.26) 



Definition 5.1. Let Cq — (Ei, E2, S3, s, cr*^, g*^, /i, iJ) be an admissible scattering con- 
figuration satisfying the following. 
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i) Si, XI2, S3 are all contained and h, H are both supported in C[^q ,p^/2,a]- 

ii) In the reg ions C j^^^^^ ^ ]\C [^Q^^^/2^a] and C[^q^^2^]\C[^q^^^^^], (a^ q^) are, respectively, 
given by (5.25) and (5.26). Moreover, in R^\C[^q ^p^,^]^ — I and — 1. 

iii) If i7 = 0, then it is required that 

^^{x) > Aq for a.e. X G supp{h), (5.27) 
where Aq is a positive constant independent of s. 

iv) If both H and /i are not identically vanishing, then it is required that 

Ao < ^q^{x),^q^{x) < Aq for a.e. x G C[^Q^^^/2,a], (5.28) 

and 

a'{x)^ ' ^ > AollCf for a.e. x G supp{H), (5.29) 
where Aq and Aq are positive constants independent of e. 
Then, we have 



Proposition 5.1. Let Cq be as in Definition 5.1, and let u^^{x,d) := Uoo{x,d;CQ) be 
the scattering amplitude corresponding to the scattering configuration Cq. Let ri ^ 1, 
/ = 1, 2 and a ^ 1 in C[^Q^a and let A/V be given in ( |4.17 ). Then there exists > 
and a function uj : (0, 60] (0, +oc] with lim5^o+ ^{^) — such that for any e < 

l|2L(-,rf)llLoo(si) <c^(^) + Ct for d^ Mr, (5.30) 

and 

\\u%o{-,d)\\L^(M^)<oj{e) + CT fordeS\ (5.31) 

where C is the constant from Proposition |4.3| and uj is independent of (Hi, 112, S3, 5). 
Moreover, ifh = and H = 0, then uj is also independent of (C[^q ,p^/2,a]5 3^^? 5^)5 
if H = 0, then uj is independent of (C[^q ,p^/2,a]5 3^^? ^5^5 ^^5^l{/i=o})5 if h ^ and 
H then uj is independent ofa^\^H=Qy 

Proof Let = (Kf , K|, , 5^ g^ /i^ i?^) be defined by 

''")*Cc- 

That is, Cq is the virtual configuration in the virtual space of the physical configuration 



Cq in the physical space. Then, by Lemma 2.1, we know Uoo{x^d;C(->) = Uoo(x^d;C(->). 



Now, the proposition can be proved by combining Theorem |3.3| and Proposition |4.3[ in 
a completely similar manner as for the proofs of Propositions |4.1| and |4.2[ □ 
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Remark 5.1. By Proposition |5.1[ we see that if r and e are chosen to be smah, the 
construction corresponding to the scattering configuration Cq yields a regularized partial 
cloak with a limited observation aperture Afr if the impinging angle is from S^, and also 
a regularized partial cloak with a limited impinging aperture Afr and if the observation is 
from S-*^. Moreover, in the cloaked region C[^q ,p^/2,a] 5 if source/sink (/i, H) is presented, 
the passive medium inside could be arbitrary (but regular); and if a source/sink term 
h is presented, one only need impose the generic condition that at the place where the 
source/sink is located, the medium is absorbing. 



Remark 5.2. The lossy layer (5.24)-(5.25) is only a very particular choice of illustrating 
our general theory. Indeed, one can devise more general lossy layers as long as As- 
sumption [!} b) is satisfied. For example, one may choose the lossy layer in (5.24) to be 



I and qf = (1 + i)e ' . Moreover, with such a lossy layer, one could set the 
parameter a in defining C[^q ^^^^j to be e. In doing this, in the virtual space, C[^q^^^^] 
is 'uniformly' small and hence the construction would yield a full regularized cloak. 
Furthermore, we would like to note that by adjusting the parameter a, one can achieve 
a customized partial cloak. That is, for any given aperture Mr and any small 5 > 0, one 



< 5 for 



can choose a sufficiently small e and a G (^,1) such that 

d G Mr- This can be straightforwardly shown by combining Theorem 3.3 and Proposi 
tion 4.3 However, there is one point that one needs to take care of in (4.20). In (4.20), 



the constant C in the estimate is dependent on a. However, heuristically speaking, C 
could be independent of a for a from a certain range. If this is the case, the construction 
of the customized partial cloak holds, and we shall not explore further in this aspect in 
the current article. 

For the three-dimensional ABC geometry D, in a similar manner, one lets 



where 



X 



Then 



D 



Identity for x E 



D 



[i(;o,^2,a] ' 



[i(;o,r2,aJ 



for X G ^\w^^e^a\ = Dl U Dm U T^r. 



[ i ^(xi + a) 

X2 
X3 




I V 



E 
£ 



\ 

/ Xl \ 



X2 



f ^ (xi — a) + a ^ 

ri 



X2 
Xs 



J 



for X G Dl, 



for X G D 



for X G Di^. 



(5.32) 



(5.33) 



'^[w,r2,a] is bi-Lipschitz and orientation-preserving and ^|aD 



[w,r2,a] 



34 



= Identity. Next, we again introduce a lossy layer as follows 

(D[^,6,a]\D[^,6/2,a]; erf , qf) , af = cie^I, ql = (c2 + ^C3)6:"^/^ (5.34) 
where Q, 1 < / < 3 are positive constants, and let 

{^[w,ri,a]\^[w,ri/2,a]'i^'' = ^* iP[w,e,a]\^[w,e/2,a]'i ^ (fl) • (5-35) 

We also let 

(D[^,^2,a]\D[^,ri,a];3^^5 5^) = ^* (D[^,^2,a] \D[^,£,a] ; l) , (5.36) 

Definition 5.2. Let C|) = (Si, S2, S3, 5, a^, g^, /i, i7) be an admissible scattering con- 
figuration satisfying the following. 

i) Si, XI2, S3 are all contained and /i, i7 are both supported in ^[w^ri/2,a]- 

ii) In the regions D[^^^^^^]\D[^^^^/2,a] and D[^^^2,a]\D[^,^,,a], (3^^5'^) are, respectively. 



given by ( |5.35| ) and ( |5.36[ ). Moreover, in R'^\D[^ ^^2,0]^ cr^ ^ I and = 1. 

iii) If i7 = 0, then it is required that 

^^{x) > Xq for a.e. X G supp{h), (5.37) 
where Aq is a positive constant independent of e. 

iv) If both H and h are not identically vanishing, then it is required that 

Ao < ^q^{x),^q^{x) < Aq for a.e. x G J^[w,ri/2,a], (5.38) 

and 

^'{x)^ • C > AollCf for a.e. x G supp{H), (5.39) 
where Aq and Aq are positive constants independent of e. 

Then, we have 



Proposition 5.2. LetCf^ be as in Definition 5.2, and let u^^(x^d) := Uooix^d^Cf^) be the 



scattering amplitude corresponding to the scattering configuration Cf^ . Letri ^ 1^ / = 1, 2 
and a ^ 1 in D[^^^ ,^^]. Then there exists > and a function uj : (0,5o] ^ (0, +oc] 
lim5^o+ ^(s) = t/iat /or any e < e^) 

Woo{-,d)\\L^{%^)<^{e) fordeS^. (5.40) 

where uj is independent of (Si, S2, S3, 5). Moreover, if h = and H = 0, then uj 
is also independent of ij^[w^ri/2,a]'i^^ 1^^)') ^'^^ if H = {), then uj is independent of 
(D[^ ,p^/2,a]5 ^5^5 ^5^l{/i=o})5 ^/ ^ ^ and H ^ 0, then uj is independent of 

^^\{H=0}' 
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Proof. The proof can be obtained in a similar manner to that for Proposition 5J_, The 
major different point is that D[^^^^^] in the virtual space degenerates to a line-segment 
as 5 ^ 0^, which has zero capacity in M^. Hence, one would have a regularized full 
cloak for the construction corresponding to Cf^; see also the proofs of Propositions |4.l| 
and 112 □ 



Remark 5.3. Similar to Remark 1 5. 2 [ the lossy layer ( 5.34 )-( 5.35) is only a very particular 



choice for illustrating our general theory. Indeed, one can also devise much more general 
lossy layers as long as Assumption [l] b) is satisfied. For example, one may choose an 
anisotropic lossy layer with variable coefficients as the one in (5.24). Moreover, we would 
like to remark that by adjusting the free parameter a to be smaller in D[^ ,p^^], one can 
improve the accuracy of approximation for the cloak construction corresponding to Cf^. 
Finally, for the ABC geometry E, one can introduce 



(Identity for x G R^\E[^Q^^2,a,6], 
^l(x) for X G ^wo,r2,aM\'^[wo,e,aM^ 

for X G E[wo,e,aM = Eo U (Uti ) , 



where 



r / 




Xi \ 

X2 



for X G Eg, 



V / 

/ ^ (xi ^ a) ± a \ 



X2 
ri 



X3 



for X e E^, 



v 

/ 

/ ^ (^1 — a) + a \ 



£ 

Xi \ 

'f{x2Tb)±b 

ri 



for X G E^, 



V 



ri 



Xs 



J 



( ^ (^1 + a) — a ^ 

^ (X2 T &) ± & 



-3^3 



for X G Eg , 



for X G Ef . 



Let a particular lossy layer be chosen to be 

[wo,e/2,a,h\ t tQi ) j 

with 



{^[wo,s,a,b]\^[wo,£/2,a,b]''' '^l ^ 1l) ' = Ci£^I, qf = (C2 + ics) 



-1/2 



(5.41) 



(5.42) 



(5.43) 
(5.44) 
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where Q, 1 < / < 3 are positive constants. Set the cloaking medium to be 



(5.45) 



Definition 5.3. Let = (Si, S2, S3, 5, a^, g^, /i, H) be an admissible scattering config- 
uration satisfying the following. 

i) Si, S2, S3 are all contained and h, H are both supported in Ej^^ ,^^/2,a,6]- 

ii) In the regions E[^o^ri,aM\^[wo,ri/2,aM and E[^Q^^2,a,6]\E[^o,ri,a,6], (3'^ q^) are, respec- 



tively, given by (5.43) and (5.45). Moreover, in 



^[wo,r2,a,b]') ^ 



I and 



iii) If i7 = 0, then it is required that 

^^{x) > Xq for a.e. x G supp{h)^ 
where Aq is a positive constant independent of e. 

iv) If both H and h are not identically vanishing, then it is required that 

Ao < "^^{x), 'R^{x) < Aq for a.e. x G 'E[^o^ri/2,aM^ 

and 

3^^(^)C • C ^ AollClP for a.e. x E supp{H), 
where Aq and Aq are positive constants independent of e. 
Then, we have 



(5.46) 



(5.47) 
(5.48) 



Proposition 5.3. Let be as in Definition 5.3, and let u^^{x^d) \— Uoo(x^d;C^) be 
the scattering amplitude corresponding to the scattering configuration C^. Let ri ^ 1, 
/ = 1, 2 and a^b ^ 1 in Ej^^^ 5] and let Mr be given in ( |4.17 ). Then there exists > 
and a function uj : (0, 60] (0, +oc] with lim5^o+ ^{^) — such that for any e < e^) 

\\u%{■,d)\\L^^s^~^<uJ{e) + CT fordeK, (5.49) 

and 

Woo{;d)\\L^^^^)<uj{£) + CT fordeS\ (5.50) 

where C is the constant from Proposition |4.3| and uo is independent of (Si, S2, S3, 5). 
Moreover, if h = and H = 0, then uj is also independent of {^[wQ^ri/2,a,h]')^^ iq^)'i 
if H = 0, then uj is independent of i^[wQ^ri/2,a,b]'i^^ ^^q^ ^^q^\{h=o})'i h ^ and 

H ^0, then uj is independent ofa^\^H=Qj. 

Proof. The proof follows in a similar manner to that for Proposition |5.1[ □ 



Remark 5.4. Similar to Remark 5.2, the lossy layer ( 5.43 )-( 5.44) is only a very particular 
choice for illustrating our general theory. One can also devise much more general lossy 
layers as long as Assumpti on [l| b) is satisfied and it could be anisotropic with variable 
coefficients as the one in (5.24). Also analogously to Remark 5.2, by adjusting the 
two parameters a and 6, one may have a customized partial cloak by the construction 
corresponding to C^. 
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6 Numerical results and discussion 



In this section, we carry out systematic experiments based on the discussions in Sec- 
tions [4] and [5] to investigate the cloaking performance with respect to the regularization 
parameter e in two and three dimensions. Numerical simulations are carried out by em- 
ploying the C++/Matlab coupled programming and the pardiso solver. Finite element 
method is used in a truncated domain coated by a perfectly matched layer to evaluate 
the performance of the proposed cloaks. The far-field pattern of the scattered field is 
calculated via the Kirchhoff-Helmholtz formula by surface integral on a sphere, within 
which the cloaking medium is compactly imbedded. 

In the following experiments, we use the L^(S^~^)-norm of the far-field patterns 
-4-T(-,rf) •= Uoo{' ^d^ C^Y^) (T = C, D, or E) to evaluate the cloaking performance of the 
specific cloaks. More precisely, we compute the module of the far field patterns on 100 
equidistant points along the unit circle in two dimensions, or 100 x 100 points on the 
unit sphere by choosing uniformly distributed polar angles and azimuth angles in three 
dimensions, respectively. The exact value of the norm of Uoo{'^d-^C^) is then approxi- 
mated by the maximum value of these modules, which is adopted for the performance 
evaluation of cloaks. 

6.1 2D cloaks 

For acoustic scattering in two dimensions, let the incident wave be a plane wave of 
wavelength 1^ — 2. We choose a — 1 ^ e to test partial cloaks and a — e to check full 
cloaks. 

First of all, we study the scattering performance of Type C cloaks. Let us visualize 
the effect of the f' and cloaks of Type C by looking at the instantaneous pressure. 
As shown by Figure [5j the total acoustic pressure is bent and/or compressed in the 
cloaking medium. There is clearly a thin surface layer within the lossy layer between 
the cloaked and cloaking regions and prevents the incident wave from penetrating into 
the cloaked region when the incident direction is parallel to the axial direction of Type 
C cloaks. At the mean time, the total pressure outside the cloaking media deviates only 
slightly from that of a perfectly plane wave (except tiny numerical deviation due to the 
finite spatial resolution in numerical simulations). 

Next, we test the performance of the 2D partial cloak in case of tilted incident 
angles. From the instantaneous pressure in Figure [6]^a), (b) and (c), it can be observed 
that there is increasingly more pressure deviation of the total field from that of a perfectly 
plane wave outside the cloak when we increase the incident angle from 5 and 10 to 20 
degrees. When the incident angle is perpendicular to the axial direction of the partial 
cloak, strong scattering amplitude occurs and this suggests the failure at the extreme 
case, see Figure [6]^d) . 

By symmetry, the same observation can be explained from the other side of the 
partial cloak. Figure [t] shows the far field pattern in decimal Bell [dB] with respect to 
the polar angle. It can be seen that when the incident direction is along the ideal axis 
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(b) (c) 

Figure 5: Real part of instantaneous acoustic pressure distribution using two- 
dimensional (a) l^, (b) P and (c) cloaks, respectively. Incident plane wave from 
left to right. Parameters: £ — 10~^, ri = 1, r2 = 2, a = 1. 



39 



(x-axis, namely or 180 degrees in polar angle), the far field is always below —40 dB[^ 
which is small and negligible in practice. Figure [6](a) suggest that for a small tilted 
incident angle below 5 degrees, approximate partial cloak effects can be achieved with 
the far-field pattern below —20 dB. 

In addition, there does exist some angle r ^ 15, 7 and 3 degrees corresponding 
to the tilted incident angles 5, 10 and 20, respectively, such that with such limited 
observation aperture perturbed off the ideal axis, the far field patterns are still below 
—40 dB (cf. cusps near the x-axis within the —40 dB contour curve in Figure [t]). Hence 
(5.31) of Proposition 5.1 is confirmed numerically. That is, the partial cloaking is still in 



effect in a limited observation aperture around the ideal direction parallel to the cloak 
even for full impinging aperture. 




I 




I 



I- 
I 



Figure 6: (From left to right) Real part of instantaneous acoustic pressure distribution 
using two-dimensional cloak with titled incident angles, 5, 10, 20 and 90 degrees, 
respectively. Parameters: £ = 10~^, ri = 1, r2 = 2, a = 1. 



^—40 dB amounts to 10 ^ in magnitude in far field patterns as opposed to the unit magnitude in the 
incident wave 
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Figure 7: Polar plot of far field pattern as a function of polar angle using two-dimensional 
cloak with tilted incident angles, 0, 5, 10 and 20 degrees, respectively. 



Finally, for the regularized full cloak case, we set a = 6: and the incident angle 
even perpendicular to the axial direction of the cloak and test the cloaking performance 
for the /2 cloak of Type C. It can be observed from Figure 8 that compared with 
the red reference line denoting the second order decay rate, the far field pattern decays 
quadratically with respect to e. That is, the construction yields an approximate full cloak 
within 6^-accuracy of the ideal full cloak. This observation complies with Remark 5.2 



Moreover, the second order rate of approximation complies with the estimate in [28] for 
regularized full cloaks. Similar observations are made for and cloaks of Type C. 



6.2 3D cloaks 

In the 3D cloaks, we fix the wavelength /g = 3 and the regularization parameter e — 10~^ 
in acoustic scattering and use the spherical domain with radius 5 centered at the origin 
as the truncated domain for numerical simulations. 

First, we check the cloaking performance of Type D cloaks in terms of a. All the 
impinging direction of the incident wave is perpendicular to the axial direction of the 
cloaks. It is noted that we always adopt the 3D cloak and the right angle incidence 
for tests. Similar observations are made for other combinations. 

We show in Figure |9](^a) the plane impinging wave with the incident direction along 
the X-axis. The resulting scattered wave is depicted in Figure[9]^b), which nearly vanishes 
outside the cloak. In Figure l9Fc), it can be obviously seen that the total wave is bent 
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Figure 8: Convergence history of scattering measurement data versus e for the 2D 
cloak when a — e. 



and compressed within the cloaking medium. 

For a full 3D cloak with a — Figure 10 demonstrates the third order approximation 
of the far-field pattern with respect to e. That is, the construction yields an approximate 
full cloak within ^^-accuracy of the ideal full cloak, which is consistent with the result in 
[28] for regularized full cloaks obtained by blowing up a 'uniformly' small region in the 
virtual space. As we increase a to be 1, the far- field patter decays only quadratically in 
terms of see Figure [TT] In other words, the construction by blowing up a 'partially' 



small region in the virtual space will suffer a reduction in the rate of approximation. 



This observation complies with Remark 5.3 



Finally, we investigate the 3D cloak of type E proposed in Definition |5.3[ We choose 
fixed length a — 1 and width h — 1 but vary the regularization parameter e. The incident 
plane wave is impinging upon the E cloak along the x-direction which is parallel to the 
xy plane where the type E cloak lies. Figure 12 shows us the scattered and total pressure 
fields in sliced plots. Obviously the scattered field is almost null outside the cloak, while 
the total field behaves nearly perfectly through the cloak by bending and compressing 
in the cloaking medium. 
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(b) (c) 

Figure 9: Type D when e — 0.01: Real part of (a) the incident plane wave expjz/cx • do}, 
(b) the scattered field and (c) the total field shced at x = 0, z = and y = 0, respectively. 
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Figure 10: Convergence history of scattering measurement data versus e for the 3D 
cloak D when a — e. 
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Figure 11: Convergence history of scattering measurement data versus e for 3D f cloak 
D when a — 1 and e varies. 
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Scattered vvave Total v^'ave 




Figure 12: Type E when e — 0.01: Real part of (a) the scattered field and (b) the total 
field sliced at x = 0, z = and y — —1.5, and 1.5, respectively. 
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